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1.  SUMMARY 


The  energy  density  in  chemical  bonds  is  a  theoretical  limit  for  most  energy  storage  systems.  For 
example,  in  fly  wheels  the  energy  density  of  the  bonds  limits  the  tensile  strength  of  the  material 
and  therefore  limits  the  energy  storage  capacity.  In  compressed  air  systems,  the  energy  density  of 
the  bonds  in  the  wall  material  limits  the  compressive  strength  and  therefore  limits  the  energy 
density  of  the  system.  The  proposed  technology  takes  advantage  of  quantization  phenomena  to 
design  nanocapacitor  storage  devices  which  reach  that  limit  and  go  beyond. 

In  Task  I  (9  months),  we  assembled  nanocapacitors  with  large  dielectric  strength  to  a  nanoscale 
rechargeable  battery.  The  following  parameters  were  investigated  experimentally:  maximum 
energy  density,  leak  currents . 

In  Task  II  (9  months)  we  fabricated  one-,  two-  and  three-dimensional  synthetic  atoms  and 
compared  their  performance  data  with  theoretical  predictions.  We  studied  recovery  after  failure 
and  likelihood  of  cascading  failures.  Simultaneously,  nanocapacitors  with  a  particularly  large 
dielectric  strength  were  fabricated  with  radio  frequency  (RF)  sputtering  in  our  laboratories.  The 
reduction  of  leak  currents  and  increased  dielectric  strength  due  to  Coulomb  blockade  was  studied 
theoretically  on  1 -dimensional  synthetic  atoms  (cathode-anode-cathode  nanocapacitors). 

In  Task  III  (9  months),  we  investigated  recovery  of  alumina  nanocapacitors  after  high-voltage 
breakdown.  Previously  we  had  discovered  that  AI2O3  nanocapacitors  can  self-repair.  Recently  we 
discovered  that  our  newly  fabricated  SiC>2  nanocapacitors  can  self-repair  too.  Nanocapacitors  with 
a  particularly  large  dielectric  strength  will  be  assembled  into  an  energy  storage  device  and  tested. 
The  data  from  task  I  were  used  to  improve  the  performance  of  the  energy  storage  device.  A  lab 
prototype  of  an  optimized  energy  storage  device  with  an  assembly  of  nanocapacitors  will  be 
completed  for  external  testing.  First  lab  prototypes  of  energy  storage  made  of  synthetic  atoms  and 
of  stacks  of  nanocapacitors  will  be  assembled  and  tested. 

In  Task  4  (9  months)  Multi-dimensional  synthetic  atoms  was  studied  theoretically. 
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2.  INTRODUCTION 


We  are  conducting  research  into  the  energy  storage  in  nanocapacitor  arrays.  Nanocapacitors  are 
expected  to  be  fully  functional  in  a  large  temperature  range,  starting  from  space  temperatures  in 
the  shadow  of  the  Earth  to  the  melting  temperature  of  tungsten  (3422  °C).  In  contrast  to  batteries 
with  electrolytes,  nanocapacitor  batteries  can  be  fully  charged  and  discharged  billions  of  times 
without  any  loss  of  efficiency.  Nanocapacitor  batteries  can  be  fully  integrated  with  advanced 
electronics.  However,  there  are  several  hurdles  to  overcome,  before  nanocapacitor  batteries  can 
replace  Li-ion  batteries  and  super-capacitors: 

(1)  The  maximum  energy  density  is  currently  much  less  than  in  Faradic  systems. 

(2)  Tunneling  currents,  field  emission  currents,  and  other  leak  currents  reduce  the  efficiency 
of  nanocapacitor  batteries. 

(3)  The  power  density  is  much  larger  than  in  Faradic  systems  and  this  is  a  potential  safety 
problem. 

We  investigate  several  hypotheses  to  overcome  these  hurdles: 

(1)  The  dielectric  strength  and  the  maximum  energy  density  depend  on  the  gap  size.  The 
dielectric  strength  of  capacitors  with  gap  sizes  in  the  nm  range  is  increased  by  3  orders  of 
magnitude.  The  energy  density  is  limited  by  the  tensile  strength  of  the  electrodes. 
Consequently,  an  array  of  a  large  number  of  nanocapacitors  has  a  much  larger  energy 
density  than  a  single,  large  capacitor  with  the  same  weight  and  volume. 

(2)  In  3 -layer,  cathode-anode-cathode  nanocapacitors  (synthetic  atoms)  leak  currents  can  be 
suppressed  with  Coulomb  blockade.  When  Coulomb  blockade  occurs,  quantum  pressure 
keeps  the  electrons  away  from  the  anode,  like  the  electron  is  kept  away  from  the  nucleus 
in  a  hydrogen  atom. 

(3)  Dielectric  nanocapacitors  can  recover  after  breakdown.  If  the  recovery  is  fast,  cascading 
failures  can  be  suppressed. 

We  are  studying: 

(1)  The  efficiency  of  energy  storage  in  nanocapacitors  arrays  as  a  function  of  temperature,  the 
number  of  charge/discharge  cycles,  and  the  average  energy  density; 

(2)  Leak  currents  in  arrays  of  synthetic  atoms  and  evidence  for  Coulomb  blockade; 

(3)  The  recovery  time  after  breakdown  and  the  probability  of  cascading  failure  in 
nanocapacitor  batteries. 

The  long-term  goal  of  this  work  is  to  create  a  new  battery  type  which  exceeds  the  energy  density, 
efficiency,  reliability  and  safety  standards  of  chemical  batteries  and  is  fully  operational  at  a  large 
temperature  range. 
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3.  TASK  I:  ASSEMBLE  NANO-CAPACITORS  WITH  LARGE 
DIELECTRIC  STRENGTH  TO  A  RECHARGEABLE  BATTERY 


3.1.  Methods,  Assumptions,  and  Procedures 

We  used  two  different  fabrication  methods.  In  one  fabrication  method  the  substrate  was  silicon.  In 
the  other  fabrication  method  the  substrate  was  glass. 

3.1.1  Fabrication  method  1 

We  bought  silicon  wafers  from  Silicon  Quest  (see  Figure  1).  The  bare  silicon  wafers  with  the 
natural  passivation  oxide  layer  are  oxidized  with  a  furnace  (see  Figure  2)  to  various  thicknesses. 


Figure  1.  Bare  silicon  wafers  with  a  natural  passivation  oxide  layer. 


Figure  2.  Photo  of  our  oxidation  furnace. 
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Then  we  checked  the  thickness  of  the  oxide  layer  with  ellipsometry  (see  Figure  3).  Finally,  we  use 
a  shadow  mask  in  the  electron  beam  evaporator  (see  Figure  4)  to  deposit  circular  metal  islands  on 
the  oxide.  Figure  5  shows  an  image  of  a  nanocapacitor,  after  using  the  electron  beam  evaporator 
to  deposit  chromium  and  gold. 


Figure  3.  This  ellipsometer  measures  the  thickness  of  the  oxide  layer. 


Figure  4.  This  electron  beam  evaporator  deposits  various  metals. 


Figure  5.  A  silicon  dioxide  layer  with  gold  island. 

The  Scanning  Electron  Microscope  (SEM)  in  Figure  6  is  used  to  determine  the  quality  of  the 
islands  (see  Figure  7). 
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Figure  6.  Hitachi  S4700  scanning  electron  microscope 


Figure  7.  SEM  image  a  nanocapacitor  (top  view  of  gold  islands 
with  diameter  1.547mm). 

We  find  that  in  the  new  samples  the  dielectric  strength  is  more  reproducible  and  larger.  Further  we 
find  that  dielectric  strength  is  larger  for  thinner  oxide  layers  (see  Figure  8).  The  dielectric  strength 
determines  the  maximum  energy  density. 
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Figure  8.  Plot  of  average  dielectric  strength  verses  thickness  of 
silicon  dioxide  nanocapacitors. 


3.1.2  Fabrication  method  2 

We  fabricated  a  multi-layer  nano-capacitor  energy  storage  systems  and  tested  them.  Figure  9  is  a 
sketch  of  the  system  on  a  glass  substrate.  The  side  view  on  the  left  and  the  top  view  is  on  the  right. 


Figure  9.  Sketch  of  3-layer  nano-capacitor  system. 
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Figure  10  shows  photos  of  the  as  fabricated  device. 


Figure  10.  Photos  of  the  3-layer  nano-capacitor  system 
described  in  Figure  1. 

A1  deposition  is  done  by  with  e-beam  evaporation.  AI2O3  deposition  is  done  in  Atomic  Layer 
Depostion  (ALD)  system.  Metal  shadow  masks  are  used  to  make  the  pattern.  The  capacitance 
was  measured  at  room  temperature.  The  amplitude  of  the  test  signal  was  0.5  V.  The  frequency 
was  120  Hz: 

Bottom-Middle  electrode  capacitance:  Cb-m=  4.62  nF 
Middle-Top  electrode  capacitance:  Cm-t=  4.18  nF 

Total  capacitance  (top-bottom  electrodes  are  connected  with  gold  wire):  Cbt-m=  8.80  nF 

The  measurements  confirmed  that  Cbt-m  =  Cb-m  +  Cm-t.  The  measured  results  are  consistent  with 
the  theoretical  estimates  if  «s4.  This  is  close  to  the  theoretical  value,  e=4.5. 

3.2  Results  and  Discussion 

We  found  that  a  stack  of  nanocapacitors  has  anticipated  energy  storage  capacity  and  extremely 
small  leak  currents.  This  might  be  due  to  Coulomb  blockade.  Preliminary  results  predict  a 
reduction  of  4  orders  of  magnitude.  The  blue  line  is  a  graph  which  includes  the  scaling  of  the 
field  enhancement  factor  with  the  thickness  of  the  dielectric.  [1] 

3.3  Conclusion 

We  conclude  that  there  are  no  apparent  obstacles  in  fabricating  a  micro-battery  with  many  layers. 
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4.  TASK  II:  FABRICATE  ARRAYS  OF  ONE-,  TWO-  AND  THREE- 
DIMENSIONAL  SYNTHETIC  ATOMS  AND  COMPARISON  OF  ITS 
PERFORMANCE  DATA  WITH  THEORETICAL  PREDICTIONS. 

4.1.  Methods,  Assumptions,  and  Procedures 

One  of  the  most  prominent  features  of  synthetic  atoms  is  charge  storage  in  the  space  between  the 
electrodes  [1,2].  Hysteresis  in  the  I-V  curve  (current  versus  voltage  plot)  is  a  footprint  of  charge 
storage  in  the  dielectric,  because  charge  originating  from  the  dielectric  has  a  larger  time  scale.  The 
I-V  curve  for  the  multi-layer  capacitor  shown  in  Figure  11  has  hysteresis.  This  supports  the 
hypothesis  that  the  theory  of  synthetic  atoms  is  a  good  model  for  such  multilayer  nano-capacitors. 


Sample  voltage  (V) 


Figure  11.  I-V  curves  of  the  3-layer  nano-capacitor  system. 


4.1.1.  Charge  Storage  between  the  electrodes 

Figure  11  shows  the  I-V  curves  of  a  3-layer  nano-capacitor  system.  The  black  and  red  curves  are 
the  I-V  curves  of  the  two  2-layer  capacitors.  The  blue  curve  is  the  I-V  curve  of  the  3-layer  system. 
The  leak  current  of  the  3-layer  system  is  not  additive  for  negative  voltages  and  particularly  small. 
These  data  seem  to  verify  the  prediction  that  synthetic  atom  store  charge  in  the  dielectric. 

For  practical  applications,  it  is  conventionally  assumed  that  most  of  the  electric  charge  in  a 
standard  metal/insulator  capacitor  is  typically  stored  on  the  metal  plates.  However,  in  a  fully 
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charged  capacitor  there  is  usually  some  extra  charge  absorbed  by  the  dielectric  media,  associated 
with  either  the  displacement  of  charges,  or  from  movement  of  electrons  from  the  metallic  plates 
of  the  capacitor  into  the  dielectric  layer  separating  the  plates.  Such  charge  penetration  into  the 
dielectric  is  driven  by  the  strong  electric  field  occurring  in  the  dielectric  spacer  of  a  nanocapacitor. 
This  absorption  of  electric  charge  by  dielectric  is  a  well-known  and  well-studied  phenomenon  [3- 
5],  Generally,  there  are  several  competing  mechanisms  by  which  the  charge  can  be  retained.  Two 
of  these  mechanisms:  hindered  dipole  rotation  and  displacement  of  ions  in  dielectric,  do  not 
involve  any  transfer  of  charge  between  the  metal  plates  and  the  dielectric  [4,  5],  Additionally,  there 
is  always  the  possibility  of  charge  penetration  into  the  dielectric  due  to  the  movement  of  hopping 
charges,  i.e.  the  tunneling  of  electrons  between  the  metal  plates  and  some  localized  charge  traps  in 
the  dielectric  [6-8]. 

Due  to  dielectric  absorption,  the  discharging  of  a  capacitor  is  essentially  a  multiple-stage  process, 
with  charge  deposited  on  metal  and  charge  acquired  by  dielectric,  occurring  on  different  time 
scales.  Obviously,  the  time  scale  for  the  discharge  of  the  plates  is  controlled  by  the  value  of  the 
shunting  resistor,  R.  Typically  the  loss  of  charge  accumulated  on  metal  plates  is  a  very  fast  process, 
described  by  the  exponential  law  with  characteristic  time  completely  determined  by  the 
capacitance  and  the  total  resistance  of  the  short  circuit.  In  contrast,  the  removal  of  charge 
accumulated  by  dielectric  is  a  considerably  slower  process,  generally  described  by  relaxation 
theory  [9], 

4.1.2.  Sample  fabrication  and  measurements 

An  aluminum  fdm  is  deposited  on  a  glass  substrate  and  covered  with  alumina  (AI2O3),  which 
serves  as  the  dielectric  layer  of  the  capacitor.  The  top  electrode  is  a  single  layer  of  graphene. 
Graphene  has  metallic  conductivity  and  acts  here  as  a  metallic  plate  of  the  capacitor.  Thus, 
effectively,  the  sample  consists  of  two  identical  capacitors,  of  the  type  Al-AbCb-graphene, 
connected  in  series.  Thus,  the  obtained  capacitance  is  twice  smaller  than  the  value  computed  for 
just  one  such  capacitor.  Note  that  the  source  of  voltage  used  here  is  adjustable,  i.e.,  it  can  be 
controlled  from  the  computer  and  can  be  set  to  zero  or  to  any  desired  value. 

The  nanocapacitors  used  in  this  study  were  made  with  a  d=50  nm  thick  A1  film  coated  by  a  7  nm 
layer  of  AI2O3  and  a  single  layer  of  graphene.  The  sample  schematic  is  shown  in  Figurel2.  A 
picture  of  the  device  is  shown  in  Figure  13.  The  length  of  the  strip  was  L=10  mm  and  the  width  is 
w=l  mm.  Thus,  the  expected  capacitance  was  C=s8oLw/d=12.6  nF.  The  actual  measured  value  was 
10  nF,  which  is  close. 
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Figure  12.  Schematic  of  a  graphene  nano-capacitor. 


Figure  13.  A  picture  of  the  device  described  in  Figure  12. 


To  measure  the  charge  Q,  stored  on  the  capacitor,  we  apply  voltage  V,  of  the  order  of  a  few  volts, 
and  wait  a  few  minutes.  Then  the  voltage  on  the  battery  is  changed  to  zero  and  the  current  is 
measured  as  a  function  of  time.  Note  that  in  equilibrium  the  current  is  zero  since  the  capacitor  does 
not  provide  a  path  for  a  dc  current.  Thus,  the  current  we  detect  is  a  transient  phenomenon,  occurring 
either  after  the  voltage  is  set  to  a  constant  nonzero  value  or  after  the  voltage  is  set  to  zero.  It  is 
interesting  to  note  that  in  a  wide  range  of  voltage  the  current  approaches  zero  with  time  even  if  the 
voltage  is  above  zero.  This  provides  strong  evidence  that  our  capacitors  are  not  leaky,  i.e.  the 
leakage  current  through  the  dielectric  is  zero. 
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4.2.  Results  and  discussion 

Figure  14  shows  I-V  curves  for  such  devices. 


U  (V)  b*U  (V) 


Figure  14.  I-V  curve  (left)  and  rescaled  I-V  curve  (right)  of  a 
graphene  nano-capacitor. 

The  capacitor  discharge  current  is  shown  in  Figure  15. 

The  initial  regime  is  the  exponential  discharge.  In  this  case  the  current  decays  exponentially  with 
time.  But  at  about  9  s  on  the  time  axis  a  new  regime  begins,  which  is  a  power-law  decay  of  the 
current.  It  is  clear  in  the  log-log  plot  shown  in  Figure  15,  where  the  second  regime  appears  to  be  a 
straight  line. 


IE-7  -a 
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Figure  15.  The  discharge  current  measured  in  the  circuit 
shown  in  Figure  12. 
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A  more  detailed  study  of  the  discharge  process  is  shown  in  Figure  16.  The  discharge  current  is 
shown  for  three  different  voltages.  In  all  three  cases  the  dependence  is  proportional  to  I(t)~l/t.  The 
exact  power  of  this  power  law  decay  can  be  either  slightly  smaller  or  slightly  larger  than  the  generic 
value  of  -l.In  this  case  the  series  resistor  was  lower  so  that  the  current  drop  occurred  faster.  So, 
when  the  measurement  was  started  only  the  power  low  residual  discharge  (stage  2)  is  observed.  It 
shows  a  power-law  dependence  quite  clearly. 


Figure  16.  The  graphene  nanocapacitor  discharge  current 

versus  time. 

It  is  interesting  that  the  integration  of  the  current,  which  gives  the  charge  stored  in  the  dielectric, 
Q  =  f  I(t)dt  gives  very  large  values,  which  are  near  the  value  of  the  charge  stored  on  the 
capacitor  plates,  QP=VC.  We  attribute  this  fact  to  the  small  thickness  of  the  dielectric,  which  is  of 
the  order  of  a  few  nanometers.  We  speculate  that  at  the  high  electric  field  created  in  the  reported 
experiments,  of  the  order  of  0. 1-1  V/nm,  the  electron  can  enter  the  dielectric  by  quantum  tunneling. 
Thus,  the  charge  penetration  is  very  massive.  For  comparison  we  note  that  in  regular  capacitors 
the  charge  stored  in  the  dielectric  is  only  a  few  percent  of  the  QP  charge. 

4.2.1.  Bulk  charge  measurements 

The  most  interesting  result  of  this  work  is  the  observation  of  a  large  charge  stored  in  the  volume 
of  the  dielectric  medium  of  the  capacitor.  We  find  that  the  charge  stored  in  the  dielectric,  Qd,  is 
approximately  equal  to  the  charge,  QP,  stored  on  the  plates  of  the  capacitor.  The  capacitor  plate 
charge  QP  is  defined  by  the  geometric  properties  of  the  plates  as  well  as  the  dielectric  constant  of 
the  dielectric,  namely  QP=VC,  where  the  geometric  capacitance  is  C=A sso/t/.  It  is  possible  to 
measure  these  two  charges,  QP  and  Qd,  separately  since  the  discharge  of  the  capacitor  is  controlled 
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by  the  resistor  R  shunting  the  capacitor  while  the  dielectric  volume  charge  Qd  decreases  at  a  much 
slower  rate,  which  is  defined  by  some  internal  process  of  the  charge  migration  in  the  dielectric. 

To  define  the  charge  Qd  we  measure  the  current  and  then  integrate  it.  The  total  charge  is  defined 
as  Q  =  Qp  +  Qd  =  f  I(t)dt  =  f  Iddt  +  f  Ipdt.  The  current  7(t)  represent  the  total  current 
flowing  from  the  capacitor.  Technically  speaking  we  cannot  measure  Id  and  IP  separately.  Yet  it  is 
clear  that  the  current  IP  can  drop  quickly,  namely  at  a  timescale  tau,  which  we  can  control  by 
choosing  the  series  resistor  R.  After  the  plates  are  discharged  almost  completely  we  can  assume 
that  /  ~  Id  and  define  the  dielectric  charge  approximately  as  Qd  =  J0  Iddt  .  Note  that  this  is  a 
lower  bound  on  the  dielectric  charge.  Yet,  as  we  see  below,  even  this  low-bound  estimate  gives  a 
dielectric  charge  which  is  approximately  equal  to  the  volume  charge. 

In  other  words,  it  is  due  to  a  discharge  of  both  QP  and  Qd.  Therefore,  to  measure  the  dielectric 
charge  we  start  the  integration  at  the  time  when  the  plate  charge  is  greatly  reduced.  The  time  scale 
of  the  discharge  of  the  capacitor  is  r  =  RC .  In  other  words,  Qp  =  CV0exp[—t/z],  where  V0  is  the 
initial  voltage.  The  corresponding  current,  as  a  function  of  time  t  is  Ip  =  (CV0/x)exp[— t/r]  = 
(V0/R)exp[-t/r]. 

4.2.2.  Discussion 

Fundamentally,  the  time-dependent  response  of  a  system  to  a  step-function  removal  of  constant 
stress  is  most  conveniently  defined  through  the  Fourier  transform  of  the  frequency-dependent 
response,  which  is  known  for  a  large  variety  of  mechanisms  of  dielectric  absorption.  As  suggested 
by  Jonscher  [9],  one  can  generally  distinguish  between  two  fundamentally  different  types  of 
dielectric  response.  The  first  type,  termed  'dipolar',  is  a  generalization  of  the  Debye  concept  of 
polarization  due  to  the  hindered  movement  of  ideal  dipoles.  This  type  of  polarization  leads  to  zero 
residual  polarization  after  discharge,  with  the  time  dependence  of  the  charge  described  by  the 
exponential  law  ~  exp(-o)pt),  where  coP  is  the  loss  peak  frequency  of  the  dielectric  response  in  the 
frequency  domain.  The  second  type  of  response  describes  the  behavior  of  a  system  where 
polarization  is  dominated  by  slow  movement  of  hopping  charge  carriers,  which  is  relevant  for  the 
charge  accumulation  due  to  tunneling  of  electrons  to  localized  states  in  dielectric.  This  type  of 
dielectric  response  has  a  fractional  power  law  time  dependence  ~  f“,  and  is  characterized  by  finite 
residual  polarization  of  the  system. 

In  practical  applications,  the  measured  property  is  usually  not  the  residual  polarization  or  the 
residual  charge,  but  the  discharge  current.  Generally,  for  times  much  longer  than  the  RC  relaxation 
time  of  the  capacitor,  the  discharge  current  due  to  different  mechanisms  of  charge  retention  by 
dielectric  can  be  described  by  the  fractional  power  law  ~  f“'7,  with  different  mechanisms  resulting 
in  different  values  of  the  power  a  [6],  The  varying  time  dependence  of  the  discharge  currents 
resulting  from  different  mechanisms  of  charge  absorption  can  in  principle  be  used  to  eliminate 
some  mechanisms  for  a  particular  dielectric,  or  to  analyze  the  relative  contributions  of  those 
mechanisms.  Additionally,  beside  the  different  time  dependences,  some  mechanisms  are 
characterized  by  dissimilarity  between  the  charging  and  discharging  processes,  which  can  be 
further  used  for  the  same  purpose.  Specifically,  the  trapped  volume  charge  formed  through  charge 
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injection,  possibly  by  tunneling  [10],  is  absorbed  and  released  on  different  time  scales,  making 
its  contribution  apparent  if  a  large  hysteresis  is  observed  between  the  charge  and  discharge 
currents. 

In  practical  terms,  the  observed  discharge  current  can  be  formally  split  into  a  contribution 
from  resistance-capacitance  (R-C)  relaxation  and  the  residual  current  due  to  dielectric  relaxation 
or  the  charge  stored  in  the  dielectric: 

I(t)  =  Io  [  e(t/RC>  ]+  i(t).  The  residual  current  can  be  further  defined  as  i(t)  =  Afa  . 

Experimentally  we  find  that  a~l.s 

4.2.3.  Model  1:  Charge  storage  in  dielectric  traps 

We  propose  a  model  which  can  reproduce  the  essential  feature  seen  in  most  experiments. 

Namely  the  model  can  predict  the  power  law  dependence  I(t)~tot’  with  a~l. 


Figure  17.  Sketch  of  the  dielectric  trap  model. 

The  model  is  illustrated  in  Figure  17.  A  plane  capacitor  is  formed  between  two  metallic  plates 
(blue)  and  a  dielectric  spacer  (yellow).  The  position  of  the  charge  traps  is  shown  by  dashed  lines. 
For  the  simplicity  of  the  argument  we  assume  that  the  charge  traps  are  positioned  along  the  dashed 
lines,  and  with  the  position  along  the  horizontal  axis  xn =dn,  where  n=l,2,3,...  is  an  integer,  and  d 
is  the  distance  between  neighbor  planes  where  the  defects  are  concentrated.  Let  us  assume  that  the 
left  metallic  plate  is  negatively  biased,  and  the  electrons  propagate  to  the  right,  by  means  of 
quantum  tunneling  [10]  (see  the  horizontal  blue  arrows).  The  tunneling  probability  to  the  first 
plane  (top  arrow)  with  localized  charge  traps  is  pl=Cexp[-kxl]=Cexp[ -kd],  the  tunneling  rate 
to  the  second  plane  (second  arrow)  is  p2 = Cexp[-kx2 ] = Cexp[-2kd /,  and  the  tunneling  rate  to 
the  third  plane  (bottom  arrow)  is  p3= Cexp[-kx3 ] = Cexp[ -3kd /,  etc.  Note  that  the  number  of 
planes  is  much  larger  than  what  we  show  on  the  schematic  of  Figure  17.  The  number  of 
planes  can  be  estimated  as  the  thickness  of  the  dielectric,  7  nm,  divided  by  the  distance 
between  the  atoms,  0.3  nm.  Such  rough  estimate  gives  about  200 planes  with  charge  traps. 

Let  us  assume  now  that  a  voltage  is  applied,  and  the  charge  begin  to  tunnel  into  the  dielectric. 
Assume  that  we  wait  time  t  and  then  estimate  how  many  planes  are  charged.  Let  us  estimate  the 
number  of  charge  trap  planes  which  can  be  charged  within  time  t.  It  can  be  said  that  the  time 
needed  to  fill  out  a  plane  n  of  charge  traps  is  Tn=M/pn,  where  M  is  the  total  number  of  charge 
traps  in  a  plane  (assumed  to  be  the  same  number  for  all  planes).  So,  a  plane  of  charge  traps 
will  be  filled  within  time  t  if  t>Tn.  The  largest  number  of  the  plane  which  will  be  filled  within 
time  t  is  Tn=t  or  pn=M/t  or  Cexp[-kxn]=M/t.  Thus  the  volume  charge  will  spread  up  to  the 
charge  trap  plane  with  the  position  xn  =-(l/k)ln[M/Ct]  or  xn=(l/k)ln[t(C/M)]=(l/k)ln[t]+C2. 
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The  total  charge  will  be  Q(t)=Mn=Mxn/d=(M/kd)ln(t)+C3 .  Finally,  the  current  is  I(t)=(M/ 
kD)r>,  as  is  observed  in  experiments. 

4.2.4.  Model  2:  Charge  storage  in  eigenstates  of  the  quantum  well  created  by  the  dielectric 

There  are  two  primary  mechanisms  by  which  charge  flows  in  a  quantum  well  under  an  electric 
field:  Fowler-Nordheim  emission  (also  called  cold  field  emission  or  CFE)  and  spontaneous 
emission.  Here,  we  consider  the  dynamics  of  the  CFE  current  of  a  discharging  nanocapacitor. 

The  cold  field  emission  current  is: 


J  =  - 


e*i  v 

0 


(7  F)2  exp 


(1) 


Where  F  is  the  macroscopic  field  strength,  &  is  the  work  function  of  the  metal,  y  is  the 
dimensionless  field  enhancement  factor,  where  is  y  «  1  for  w  <  30  nm.  v  «  1  the  dimensionless 
correction  coefficients  determined  from  a  numerical  analysis  of  the  surface  of  the  metal  (or  a  more 
careful  analytical  treatment),  ai  and  ai  are  the  Fowler-Nordheim  constants. 


Consider  a  two-plate  capacitor  separated  by  a  width  of  w  and  with  a  potential  difference  of  A  V. 
The  field  strength  between  the  plates  is  F  =  A  V/w  and  the  charge  density  on  the  plates  is  a  =  e  F. 
If  the  capacitor  is  disconnected  from  the  potential  difference,  current  J  will  flow  which  will  drive 
a  change  in  the  field  strength.  Thus,  the  previous  equations  can  be  rewritten  as  a  differential 
equation  for  the  field  strength  noting  that  for  two  parallel  plates:  J(t)  =  eF’(t) 


F\t)  =  (7F(t))2exp(^— 

e0  V  7 F(t)  J 


(2) 


This  Ordinary  Differential  Equation  (ODE)  can  be  solved  analytically  with  the  initial 
condition  F(0)  =  Fo  and  solutions  take  the  form: 
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<p3/2\ 


I  log 


exp 


«2 


03/2 ' 


jFo 


+ 


Qq  Ct!2  7 


01/2 


,  -1 


(3) 


We  find: 
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This  expression  for  the  current  (roughly)  exhibits  the  J~  1  / 1  behavior  as  previously  measured. 
In  the  other  case,  where: 


( «2  '#>3/2  N\  _  _  «l  «2  7  ^  >A1/2  . 

oxplTfTj  ^ - « - f 

The  current  is  approximately  constant: 

T t  j.\  _  "1 72  Fq  v  (  «2  4>Z/2  \ 

J{t)  *  — — exp 


(5) 


(6) 


Figure  18  shows  a  typical  dynamics  of  the  cold  field  emission  current  versus  time  for  an  initial 
field  strength  of  Fo  =  1  eV/nm. 

Thus,  there  is  a  "critical"  time,  rc,  when  the  current  transitions  from  being  relatively  constant  to 
falling  off  like  1  / 1: 


rc 


6 

c*i  c*2  7 v  (ft1/2 


exp 


/  a2  4>Z/2  \ 
V  7  Fo  ) 
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CFE  Current  Dynamics 


t>Tc 

t  >  tc  1  st  order 
t<  Tc 


Figure  18.  The  CFE  current  as  a  function  of  time  and  the 
approximation  functions. 

Figure  19  shows  the  break  down  time  versus  the  field  enhancement  facto  for  typical  parameter 
values. 
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Field  Enhancement  Factor  vs.  Breakdown  Time 


Y 

Figure  19.  The  breakdown  time,  xc  versus  field  enhancement 

factor  Y- 


4.2.5.  Estimation  of  the  leakage  current 

Based  on  the  model  presented  above,  and  the  fact  that  the  power  law  behavior  of  the  current  was 
observed  at  all  time  scales  tested,  one  should  assume  that  the  electrons  cannot  tunnel,  within  the 
time  scale  of  the  experiments,  through  the  entire  thickness  of  the  dielectric.  Thus,  the  leak  current 
should  be  very  small.  In  order  to  determine  the  leakage  current  experimentally  we  have  charged 
our  capacitors  up  to  3  V,  then  disconnected  the  capacitor  from  the  charging  circuit  and  allowed  it 
to  sit  without  being  connected  to  anything  for  time  At.  In  this  time  interval  the  discharge  of  the 
capacitor  was  happening  only  through  the  leakage  current  in  the  dielectric.  The  results  are  shown 
in  Figure  20.  The  discharge  is  done  not  immediately  after  the  capacitor  was  charged  but  after  some 
delay  t2.  The  delay  for  the  plotted  curves,  from  the  top  to  the  bottom,  was  t2=2  min,  10  min,  1.5  h, 
4h,  23  h.  The  exponential  discharge  fit  is  done  for  the  case  of  t2=4h.  The  initial  current  was  22.3 
nA,  according  to  the  fit. 
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Figure  20.  The  discharge  current  of  a  capacitor  through  a 
resistor  of  Rs=100  MOhm. 

The  initial  charge  was  Qi=CVi.  The  initial  voltage  is  Vi=Ii*Rs,  where  the  shunting  resistor  is 
Rs=100  MOhm.  The  initial  current  is  22.3  nA.  so,  the  initial  voltage  is  Vi=2.23  V.  Thus  within  4 
hours  the  charge  dropped  by  AQ=C(Vo-Vi)  =8.4  nC.  So,  the  leak  current  roughly  estimated  is 
I_leak=AQ/At=0.6  pA.  With  this  level  of  leakage,  the  discharge  time  is  a  few  hours.  Further 
improvements  of  the  dielectric  will  be  needed  to  store  charge  over  longer  time  intervals.  In  the 
future  we  plan  to  test  the  effect  of  low  temperatures.  It  is  expected  that  leakages  will  be  reduced  at 
low  temperatures. 

4.2.6.  Total  charge  stored  in  the  dielectric 

Here  we  present  an  example  of  our  data  analysis  intended  to  define  the  amount  of  charge  stored 
on  the  dielectric.  The  charge  on  metal  (total  charge  from  the  exponential  current)  is  well  defined, 
simply  because  it  can  be  fitted  with  an  exponential  function  and  then  extrapolated  to  zero  time  and 
to  the  infinity.  Such  extrapolated  exponent  can  then  be  easily  integrated  mathematically  since  the 
integral  of  the  exponent  is  finite.  On  the  other  hand,  there  is  much  noise  in  the  measured  current 
data  so  that  the  overall  charge  is  less  precisely  defined.  Although  the  noise  is  integrated  out  to  zero, 
since  it  is  Gaussian  and  the  fluctuations  occur  to  the  positive  direction  with  the  same  probability 
as  the  negative  direction.  On  the  other  hand,  there  is  another  inaccuracy,  related  to  the  fact  that  the 
measurement  apparatus  (Keithley  65 1 7)  gives  some  offset.  To  eliminate  the  offset  contribution  we 
determine  the  current  at  very  large  times.  Then  we  subtract  this  offset  current  from  the  raw  data, 
only  after  that  the  integration  is  done.  Thus,  the  obtained  estimate  to  the  stored  charge  is  Qd=3. 0 
10'8  C=30  nC.  On  the  other  hand,  the  charge  stored  on  the  metallic  plates  of  the  capacitor  is  QP  = 
CV.  The  applied  voltage  was  V=3  V.  The  capacitance  was  estimated  from  the  discharge  curves  and 
was  found  C=12.5  nF.  So  the  charge  on  the  plates  is  QP=37.5  nC.  We  find  that  Qd  is  just  slightly 
lower  than  QP.  Thus,  the  results  indicate  that  the  dielectric  can  double  the  efficiency  of  the  energy 
storage  in  nano-capacitors.  The  total  stored  charge  is  67.5  nC  according  to  our  analysis.  Thus,  the 
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total  energy  is  E=QPV/2=100  nJ.  The  energy  density  is  E/A/dl=48  MJ/m3.  This  can  be  converted 
to  the  energy  per  unit  mass.  The  density  of  AI2O3  is  3.95  g/cm3,  so  the  density  is  3950  kg/m3. 
Thus,  the  mass  of  the  dielectric  is  8.4  nano-gramm.  So,  the  energy  density  is  12  MJ/kg.  This  is 
about  four  times  less  than  for  gasoline.  Still,  this  result,  taken  together  with  the  observed  slow 
discharge,  suggests  that  nano-capacitors  can  be  used  for  energy  storage. 

4.2.7.  Three-dimensional  synthetic  atoms:  The  point  of  closest  distance  of  two 
perpendicular  plates  which  almost  touch 

Two  perpendicular  gold  foils  are  attached  to  a  commercially  available  nano-positioning  device. 
The  device  pushes  the  foils  towards  each  other.  When  an  electrical  contact  is  established,  the 
motion  is  reversed,  to  establish  a  gap  of  a  desired  size.  A  typical  gap  is  of  the  order  of  -100  nm, 
so  it  is  not  visible  on  the  photograph. 

We  measured  break  down  voltage  as  a  function  of  the  distance  between  gold  plates  (see  Figure  21 
and  Figure  22).  The  edges  were  electropolished  after  mechanical  polishing  before  measurements. 
Applied  voltage  was  removed  immediately  after  breakdown.  The  approach  of  edges  was  stopped 
as  soon  as  they  touched  each  other.  Then  the  plates  were  pushed  to  the  desired  distance. 


Gold  Foil 


Figure  21.  Schematics  of  a  3D  synthetic  atom:  Two 
perpendicular  gold  foils. 

The  experiments  (see  Figure  22)  suggest  that  vacuum  capacitors  can  reach  an  average  electric  field 
of  2.5  GV/m,  which  is  a  factor  of  5  larger  than  dielectric  capacitors.  We  assume  that  the  dielectric 
material  has  a  larger  tensile  strength. 
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Figure  22.  The  breakdown  voltage  (left)  and  the  corresponding 
dielectric  strength  versus  the  distance  (right). 

4.2.8  Two-dimensional  synthetic  atoms:  A  Spindt  ridge 

Next  we  designed  a  two-dimensional  synthetic  atom  with  a  linear  cathode  and  a  planar  cathode 
(see  Figure  23-25).  The  anode  is  a  Spindt  ridge  in  a  trench  in  a  SiCh  layer  and  the  counter  electrode 
is  a  planar  graphene  electrode  (Wire  width  -130  nm  Spindt  tip  height  -200  nm  Distance  from  the 
tip  to  the  graphene  -100  nm). 


Graphene 


Si3N4  (60  nm) 

SiOz  (500  nm)  ^ 
Si  (0.5  mm) 

©  °  © 

6  ^ 

'  Spindt  tip 


Figure  23.  Schemtics  of  a  2D  synthetic  atom:  A  Spindt  ridge  in 

a  trench 
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Figure  24.  SEM  images  of  a  Spindt  ridge  in  a  trench  (left)  and 

the  top  view. 


Figure  25.  Picture  of  a  Spindt  ridge  in  a  trench 

Figure  26  is  the  top  view  the  trench  without  and  with  the  Gold  Spindt  ridge.  The  gold  on  the  flat 
surface  is  eventually  pulled  off. 

Figure  27  shows  the  device  with  a  graphene  layer. 

Figure  28  shows  typical  I-V  curves  of  the  graphene  nano  capacitor  at  room  temperature.  The  leak 
current  is  much  smaller  when  the  charge  of  the  Spindt  ridge  is  positive. 
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Figure  29  shows  typical  I-V  curves  of  the  graphene  nano  capacitor  at  77K.  The  I-V  curve  is  similar 
to  the  I-V  curve  at  room  temperature,  but  has  more  features.  We  are  working  on  a  theory. 
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Figure  26.  Trench  before  and  after  Au  deposition. 
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Figure  27.  Sample  images  after  graphene  deposition 
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Figure  28.  Current- Voltage  characteristics  at  room 
temperature  in  air. 


Figure  29.  Current- Voltage  characteristics  at  77K. 


Figure  30  shows  damage  done  to  the  device  due  to  break  down.  Figure  31  is  a  device  with  a 
discontinuous  Spindt  ridge  due  to  a  stitching  problem  with  e-beam  lithography. 
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Figure  30.  SEM  image  of  the  sample  after  local  breakdown. 


Figure  31.  Bridges  in  trenches  create  interuptionds  of  the 

Spindt  ridge. 

In  summary,  the  I-V  curves  (Current  versus  Voltage  curves)  at  room  temperature  match 
expectation,  but  are  more  complex  than  expected  at  77K. 

Figure  32  shows  attempts  to  control  bridges  in  the  trenches  and  other  imperfections. 
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(d) 


(e) 


Figure  32.  Controlling  bridges  and  stiching  problems  with 
angled  arrays  of  trenches  (a-e) 


4.2.9.  One-dimensional  synthetic  atom:  A  parallel  plate  vacuum  nano-capacitor 

We  started  to  fabricate  parallel  plate  vacuum  nano-capacitors.  Figure  33  shows  the  top  view  of 
nano-vacuum  capacitors  with  Si02  spacers. 
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Figure  33.  Fabrication  of  arrays  one-dimensional  synthetic 

atoms. 
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5.  TASK  III:  INVESTIGATE  RECOVERY  OF  ALUMINA 

NANOCAPACITORS  AFTER  HIGH-VOLTAGE  BREAKDOWN 

5.1.  Methods,  Assumptions,  and  Procedures 

Previously  we  had  discovered  that  AI2O3  nanocapacitors  can  self-repair  [1],  Recently  we 
discovered  that  our  newly  fabricated  SiCh  nanocapacitors  can  self-repair  too.  Figure  34  shows  in 
red  the  applied  voltage  across  a  sample  and  a  10  MegaOhm  resister  in  series.  The  blue  curve  is  the 
sample  resistance.  This  figure  illustrates  that  SiCh  nano-capacitors  recover  after  breakdown. 

Figure  34  shows  the  resistance  versus  time  of  a  113.5nm-thick  Si02  nanocapacitor.  This  plot  shows 
a  breakdown  on  the  first  increase  of  voltage  and  shows  that  when  the  voltage  was  applied  at  0.2V 
or  -0.2V  for  an  hour,  the  resistance  increases.  Since  the  resistance  increases  at  low  voltages  after 
dielectric  breakdown  we  can  conclude  that  the  dielectric  recovers. 


Figure  34.  Plot  of  resistance  on  sample  and  applied  voltage 

verses  time 
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5.2.  Results  and  Discussion 

The  recovery  is  on  the  minute-timescale.  This  suggests  that  it  is  a  chemical  process.  We  think  that 
at  breakdown  the  Aluminum  electrode  forms  a  thin  wire  from  the  cathode  to  the  anode.  If  the  wire 
has  a  diameter  of  less  than  8  nm  it  will  fully  oxidize  in  the  oxide  rich  environment  of  the  dielectric 
layer  within  a  few  minutes.  Therefore,  the  aluminum  layer  has  to  be  in  the  10  nm  range.  Otherwise 
the  wire  is  thicker  and  does  not  fully  oxidize  after  breakdown. 

5.3.  Conclusion 

This  means  that  nanocapacitors  are  self-repairing  after  breakdown,  if  the  aluminum  electrode  has 
a  thickness  in  the  10-nm  range  and  the  dielectric  is  in  the  10-nm  range. 

We  submitted  a  paper  entitled  “Recovery  of  alumina  nanocapacitor  after  high-voltage 
breakdown”  [11],  We  filed  a  disclosure. 
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6.  TASK  IV:  INVESTIGATE  THE  REDUCTION  OF  LEAK  CURRENTS 
AND  INCREASED  DIELECTRIC  STRENGTH  DUE  TO  COULOMB 
BLOCKADE  WILL  BE  STUDIED  THEORETICALLY  ON  MULTI¬ 
DIMENSIONAL  SYNTHETIC  ATOMS 

A  3-dimensional  synthetic  atom  is  a  spherical  capacitor  on  a  nano-scale. 

The  electric  field  for  a  spherical  capacitor  is  given  by: 

Q 

E(r)  - - rfor/?i  <  r  <  R2. 

UTfoTp-H  ^ 

To  solve  the  Schroedinger  equation  for  a  spherical  capacitor,  we  make  the  following  assumptions: 
We  assume  that  there  is  only  one  electron  in  the  capacitor  that  is  we  assume  that  there  is  no 
interaction  between  the  electrons. 


h2  „  \ 

-—V2  +  V(r)  'I'  =  EH. 

2m  / 

For  a  spherical  capacitor  the  potential  is  given  by 

r  <  Ri 
Ri  <  r  <  R2 
R-2<  r  <  R:i 
R,<r 


(9) 


(10) 


Where  Ki  and  K2  are  constants  (related  to  the  Fermi  level)  dependent  on  the  material  used.  For 
now,  we  assume  that  the  potential  goes  to  zero,  if  r  goes  to  infinity,  where  R3  is  the  radius  of  the 
entire  capacitor,  including  the  outer  shell.  One  can  show  that  the  momentum  operator  is 
independent  of  r.  There  on  can  separate  the  wave  function  into  a  spherical  part  and  a  radial  part 
y/=  R(r)Y(0,  <P).  So  we  get  for  the  Schroedinger  equation: 

-  H~-d?{rRY)  +  t^t(  W)  +  V(t)RY  =  ERY 

2m  r  2mr~ 


We  can  separate  the  spherical  part  Y(0,  0)  and  u(r)=r  R(r) 


fi¬ 

lm 


t/(r)  +  ^//(rMr)  -  Eu(r). 


(12) 
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Where  we  use 


V.„  =  V(r)  +  (] 

Results  with  the  values  given:  E=  1  1  V/nm  and  r=lnm.  We  can  use  the  formula  for  the  electric 
field  we  get:  E  =Q/  (4neor2).  Solving  this  equation  for  Q  we  get  Q=1.11107xl0-19C.  For  the 
distane  between  the  first  two  energynlevels  n=l  and  n=2: 


El-Ei  = 


(1  -  1/4)  =  —  4.9ileV 


(14) 


For  the  value  of  kBT  at  room  temperature  T=298K  we  get 

knT  =  2.5710-2eV  (15) 

So,  we  see  that  quantization  happens  at  room  temperature  because  the  value  of  IcbT  is  a  hundred 
times  smaller  than  the  distance  between  the  first  two  energy  levels. 

The  wave  function  for  a  hydrogen  atom  (which  is  the  same  for  our  simplified  model  if  the  charge 
in  the  center  is  Q  =  e)  is  given  by  [1]: 


tn-(-l)'. 

2 n(n  -F  £)! 


(16) 


p  is  given  by  p=2r  /  (n  ao).  If  we  want  to  adjust  this  wave  function  for  a  center  charge  Q,  we 
can  use  the  definition  of  the  Bohr  raduis  ao  =  ( 4neoh 2)  /  (mee2)  and  exchange  one  electron  with  the 
charge  on  the  inner  sphere.  If  a  o  denotes  the  “Bohr  radius”  of  the  capacitor,  we  can  write  ao  =  ao 
Q/e  .  This  way  we  can  easily  adjust  the  wave  function  and  any  other  problem  for  our  case. 


Now  let  us  take  a  look  at  the  degeneracy  of  the  states.  The  quantum  number  1  is  defined  as  0  <  1  < 
n-1  and  the  quantum  number  is  defined  by  0  <  |m|  <  1.  Therefore  we  have  for  each  1-value 
21+1  degenerate  states  and  for  each  n  we  have: 


r=o 


(17) 


degenerate  states.  If  we  assume  that  the  wave  function  behaves  as  the  wave  funtion  for  a 
hydrogen,  we  can  calculate  the  number  of  electrons  if  all  states  are  occupied.  We  can  put  two 
electrons  in  each  state  due  to  Pauli’s  exclusion  principle,  For  the  radius  of  one  orbital  we  use  the 
following  formula: 
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So  we  get  for  the  first  orbital  inside  the  capacitor  r=Ri  and  the  last  one  that  is  inside  the 
capacitor  r  =  Rn  : 


f,  7°^ 

Rl  "■‘IT 


(19) 


So  we  get  for  the  number  of  electrons  in  the  capacitor: 


n*  ” 2  x  it  ”3 = - n?) + 

n*ni 


+  +  ^{"2  -  «|J 


(20) 


We  know  that  m  is  proportional  to  (Ri)1/2  (we  assume  that  we  get  integer  values),  we  can  now 
compare  this  quantity  to  the  volume  of  the  capacitor  V  =  4/3  n(R23  -Ri3). 


We  can  calculate  the  transition  rate  between  two  states  using: 


Sircofic3  (21) 

Q  denotes  the  emission  frequency,  d=e2|<l|r|2>|2  is  the  mean  squared  average  of  the  electric 
dipole  moment  [2][3].  The  emission  frequency  is  given  by  o=  En2  -  Eni.  Using  these  results  we 
can  now,  for  example,  calculate  the  transition  rate  for  2P  <-2S.  To  get  the  dipole  moment,  we 
implemented  the  wave  functions  in  Mathematica.  Additionally  we  compted  the  trasition  rates  for 

<l|n|2>. 

We  calculated  the  transition  rate  for  5P  ->4S,  because  the  n=4  orbital  is  the  first  that  is  actually 
inside  the  capacitor.  So  this  is  between  the  two  lowest  states  of  the  capacitor.  The  transition 
rate  is  independent  of  m,  we  used  m=0.  For  the  calculations  we  used  for  the  charge  of  the  inner 
sphere  Q  =  e,  but  the  calculations  could  be  easily  redone  for  a  different  charge: 


(1.  (I.  0|j|5,  1,0)  =  0 
(4, 0. 0|jr/|5, 1,0)  =0 
{4,0,01*15,1,0}  =  4,918  do 
d 2  —  e2al  x  24.185 


(22) 
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And  we  get  for  the  emission  frequency: 


E\  =  — O.SSeV 
£s  =  — 0.544eV 

AE  =  Eh  -  Ea  =  hu  <=>  ta  =  4.65  x  10u  s_t. 


And  finally  get  for  the  transition  rate: 


Ar. 


’ iP-t  IS 


a-E0ht'1 


=  7.37  x  lO8*-1. 


The  quantum  capacitance  for  a  spherical  quantum  dot  is  given  by: 


(23) 


(24) 


Cq(N)  = 


E( N  +  1)  -  2 E(N)  +  E(X  -  1)‘ 


(25) 


The  quantum  capacitance  depends  on  the  occupied  states.  N  is  is  the  Nth  particle  of  the  system. 
E(N)  is  the  energy  of  that  particle  [6],  The  theoretical  model  that  describes  quantum  is  given  below. 

[7]: 


/i(.V)  —  E{N)  —  E(N  —  1). 


(26) 


The  chemical  potential  describes  the  energy,  which  we  need  to  add  to  one  electron. 


1  _  AV 
C~AQ 

AV,  =  ^  _  ).(-V  +  A.V)  -  ,.(.Y) 
e 

We  set  AN  —  1 


=>Cq(N) 


E(N  +  1)  —  2E(N)  +  E(N  —  1)' 


(27) 


These  calculation  suggest  that  the  leakage  current  at  room  temperature  is  very  small  in  a 
spherical  nanocapacitor  due  to  Coulomb  blockade. 
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7.  CONCLUSIONS 


From  a  research  perspective,  the  largest  achievement  of  this  project  was  that  Coulomb  blockade 
describes  the  behavior  of  the  charge,  i.e.  there  are  space  charges  inside  the  dielectric  which  need 
to  tunnel  through  a  surface  barrier  and  surface  charges  which  leave  the  surface  on  a  time  scale 
given  by  the  external  resistor  and  the  capacity  of  the  capacitor.  In  our  nanocapacitors,  the  time 
scale  of  the  space  charges  is  much  larger  than  the  time  scale  of  the  surface  charges. 

From  a  practical  perspective,  the  most  important  finding  is  that  there  seem  to  be  no 
obstacles  in  fabricating  stacks  of  nanocapacitors  with  many  layers. 
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LIST  OF  SYMBOLS,  ABBREVIATIONS,  AND  ACRONYMS 


A1 

Aluminum 

AI2O3 

Aluminium  oxide,  Alumina 

ALD  system 

Atomic  Layer  Depositon  system 

CFE 

Cold  Field  Emission 

I-V  curve 

A  graph  of  the  current  I  versus  the  voltage  V 

ODE 

Ordinary  Differential  Equation 

RC 

Resistance-Capacitance 

RF 

Radio  Frequency 

SEM  image 

Scanning  Electron  Microscope 

Si02 

Silicon  dioxide 
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